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The invariant differential cross section for inclusive electron production in p + p collisions at 
\fs = 200 GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion 
Collider over the transverse momentum range 0.4 < pr < 5.0 GeV/c at midrapidity (\r]\ < 0.35). 
The contribution to the inclusive electron spectrum from semileptonic decays of hadrons carrying 
heavy flavor, i.e. charm quarks or, at high pr, bottom quarks, is determined via three independent 
methods. The resulting electron spectrum from heavy flavor decays is compared to recent leading 
and next-to-leading order perturbative QCD calculations. The total cross section of charm quark- 
antiquark pair production is determined as a C c = 0.92 ± 0.15(stat.) ± 0.54(sys.) mb. 

PACS numbers: 13.85.Qk, 13.20.Fc, 13.20.He, 25.75.Dw 



The production of hadrons carrying heavy quarks, i.e. charm or bottom, serves as a crucial proving ground 
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for quantum chromodynamics (QCD), the theory of the 
strong interaction. Because of the large quark masses, 
charm and bottom production can be treated by pertur- 
bative QCD (pQCD) even at small momenta without be- 
ing significantly affected by additional soft processes 0. 
This is in distinct contrast to the production of particles 
composed solely of light quarks, which can be evaluated 
perturbatively only for sufficiently large momenta. Con- 
sequently, pQCD calculations of heavy quark production 
are expected to be reliable over the full momentum range 
experimentally accessible at collider energies. 

For bottom production, next-to-leading order (NLO) 
calculations arc in reasonable agreement with data Q. 
Charm measurements at ^/s = 1.96 TeV exist for high 
transverse momentum (pr) only where the cross sec- 
tion is higher than NLO predictions by > 50%. However, 
these discrepancies are within the substantial experimen- 
tal and theoretical uncertainties 0. At the Relativistic 
Heavy Ion Collider (RHIC), charm data have been shown 
for p+p and d + Au collisions at ^/snn = 200 GeV 
as well as for Au+Au collisions at 130 and 200 GeV f 
Further measurements are crucial for a better under- 
standing of heavy flavor production at RHIC. In par- 
ticular, the relevance of higher order processes and other 
production mechanisms like jet fragmentation is unclear. 

We report on the midrapidity production (\r]\ < 0.35) 
of inclusive electrons, (e + + e _ )/2, in p + p collisions 
at y/s = 200 GeV measured by the PHENIX experi- 
ment |8( at RHIC. Contributions from semileptonic heavy 
flavor decays are extracted in the electron pr range 
0.4 < px < 5.0 GeV/c. The resulting invariant differen- 
tial cross section is an important benchmark for pQCD 
calculations of heavy quark production. Furthermore, it 
provides a crucial baseline for measurements in nuclear 
collisions at RHIC. Since hadronic heavy flavor produc- 
tion is expected to be dominated by initial parton scat- 
tering, systematic studies in p + p and d + Au collisions 
should be sensitive to the nucleon parton distribution 
functions as well as to nuclear modifications of these such 
as shadowing Q. In Au+Au collisions, heavy quarks con- 
stitute a unique and, with the data presented here, cali- 
brated probe for the created hot and dense medium. Pos- 
sible medium effects on heavy flavor probes include en- 
ergy loss 0, ^3 , az imuthal asymmet ry JT3l . and quarko- 
nia suppression [p| or enhancement [lj, [l5| . 

The data used here were recorded by PHENIX during 
RHIC Run-2. Beam-beam counters (BBC), positioned at 
pseudorapidities 3.1 < \r)\ < 3.9, measured the collision 
vertex and provided the minimum bias (MB) interaction 
trigger defined by at least one hit on each side of the ver- 
tex. Events containing high pr electrons were selected by 
an additional level- 1 trigger in coincidence with the MB 
trigger. This lcvel-1 trigger required a minimum energy 
deposit of 0.75 GeV in a 2 x 2 tile of towers in the elec- 
tromagnetic calorimeter (EMC) 0|. After a vertex cut 
of \z v tx\ < 20 cm, an equivalent of 465 x 10 6 MB events 



sampled by the EMC trigger was analyzed in addition to 
the 15 x 10 6 events recorded with the MB trigger itself. 

The PHENIX east arm spectrometer (\r)\ < 0.35, 
A<fi = 7r/2) includes a drift chamber and a pad cham- 
ber layer for charged particle tracking. Tracks were con- 
firmed by hits in the EMC matching in position with the 
track projection within 3cr. Electron candidates required 
at least two associated hits in the ring imaging Cerenkov 
detector (RICH) in the projected ring area. Random co- 
incidences of hadron tracks and hits in the RICH occured 
with a probability of (3.0± 1.5) x 10 -4 . For electrons the 
energy E deposited in the EMC equals the momentum p. 
Requiring \(E — p)/p\ < 3cr, a total charged hadron rejec- 
tion factor of about 10 4 (10 5 ) was achieved for pr = 0.4 
(> 2.0) GeV/c. Remaining background (< 1 %) was 
measured via event mixing and subtracted statistically. 

The differential cross section for electron production 
was calculated as 
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where J Cdt is the integrated luminosity measured with 
the MB trigger or sampled with the EMC trigger, re- 
spectively, tMas is the probability for an electron event 
to fulfill the MB trigger condition, N e is the measured 
electron yield, and Ae rec is the product of geometrical 
acceptance and reconstruction efficiency. For the EMC 
triggered sample, e rec includes the trigger efficiency e\ v i\. 

J Cdt is calculated as Nmb/pbbc, where Nmb is the 
number of MB triggers or, for the EMC triggered sample, 
the number of EMC triggers divided by the measured 
fraction of MB events which simultaneously fulfill the 
EMC trigger criterion. With the MB trigger cross section 
&bbc = 21.8 ± 2.1 mb 16], the analyzed data samples 
correspond to integrated luminosities of 0.7 nb _1 (MB 
trigger) and 21 nb _1 (EMC trigger), respectively. Thepr 
independent trigger bias e^ias = 0.75±0.02 was measured 
for events containing a tt withp T > 1.5 GeV/c (T^ and 
confirmed for charged hadrons with pr > 0.2 GeV/c [l7j . 
indicating a universal bias both for hard and soft pro- 
cesses. Ae rec was calculated as a function of pr (< 10 % 
variation over the full pr range) in a GEANT |l8| sim- 
ulation of electrons with flat distributions in rapidity 
(\y\ < 0.6), azimuth (0 < <j) < 2tt), and event vertex 
(\z\ < 30 cm) as input. The simulated detector response 
was carefully tuned to match the real detector. Rigor- 
ous fiducial cuts were applied to eliminate active area 
mismatches between data and simulation as well as run- 
by-run variations. The trigger efficiency ti v n, evaluated 
for single electrons in the fiducial area, rises from zero at 
low p T to 95 ± 5% for p T > 2 GeV/c. Finally the effect 
of finite bin width in pr was appropriately corrected for. 

The corrected electron spectra from the MB and EMC 
triggered samples cover pr ranges of 0.4 < pr < 
2.0 GeV/c and 0.6 < p T < 5.0 GeV/c, respectively. They 
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FIG. 1: (Color online) Inclusive electron invariant difTerential 
cross section, measured in p + p collisions at yfs = 200 GeV, 
compared with all contributions from electron sources in- 
cluded in the background cocktail (upper panel). Error bars 
(boxes) correspond to statistical (systematic) uncertainties. 
Relative contributions of all electron sources to the back- 
ground cocktail (lower panel). 



are consistent with each other within the statistical un- 
certainties in the pt region of overlap. The weighted 
average of both measurements is shown in Fig. ^ 

The systematic uncertainty of the inclusive electron 
spectrum is about 12%, almost pt independent, calcu- 
lated as the sum in quadrature of contributions from the 
acceptance calculation (7%), electron identification cuts 
(5.2%), run-by-run variations (4%), tracking efficiency 
(3%), momentum scale (1 - 5%), and other smaller un- 
certainties. The value of 12% does not include the 9.6% 
uncertainty of the absolute normalization. 

The invariant cross section of electrons from heavy fla- 
vor decays was determined by subtracting a cocktail of 
contributions from other sources from the inclusive data. 
The most important background is the 7r° Dalitz decay 
which was calculated with a hadron decay generator us- 
ing a parameterization of measured tt° Q and *± 
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FIG. 2: Ratio of electrons from heavy flavor decays (non- 
photonic) and other sources (photonic) , Rnp, for three inde- 
pendent analysis methods. Error bars (boxes) are statistical 
(cocktail systematic) uncertainties. 



spectra as input. The spectral shapes of other light 
hadrons h were obtained from the pion spectra by mr 
scaling. Within this approach the ratios H/tc are con- 
stant at high pt and for the relative normalization we 
used: r//ir a = 0.45 ± 0.10 [H, p/tt° = 1.0 ± 0.3, uj/it = 
1.0 ± 0.3, V/tt = 0.25 ± 0.08, and = 0.40 ± 0.12. 

Only the r\ contribution is of any practical relevance. An- 
other major electron source is the conversion of photons, 
mainly from 7r° — > 77 decays, in material in the accep- 
tance. The spectra of electrons from conversions and 
Dalitz decays are very similar. In a GEANT simulation of 
7T° decays, the ratio of electrons from conversions to elec- 
trons from Dalitz decays was determined as 0.73 ± 0.07, 
essentially px independent. Contributions from photon 
conversions from other sources were taken into account 
as well. In addition, electrons from kaon decays (K e3 ), 
determined in a GEANT simulation based on measured 
kaon spectra [l9j |. and electrons from external as well 
as internal conversions of direct photons 0, H3| were 
considered in the cocktail. All background sources are 
compared with the inclusive data in the upper panel of 
Fig. n with the relative contributions shown in the lower 
panel. The total systematic uncertainty of the cocktail 
is about 12%, essentially pr independent. This uncer- 
tainty is dominated by the systematic error of the pion 
parameterization (« 10%). Other systematic uncertain- 
ties, mainly the 17/ 7r° normalization and, at high px, the 
contribution from direct radiation, are much smaller. 

Given the small amount of material in the acceptance 
(Be beam pipe: 0.29 % A ; air: 0.28 % X ) the ratio 
Rnp of non-photonic electrons from heavy flavor decays 
to background from photonic sources is large (Rnp > 1 
for px > 1.5 GeV/c) as shown in Fig. [3 Two comple- 



5 



mentary analysis methods confirm the cocktail result: 

The converter technique compares electron spec- 
tra measured with an additional photon converter Xc = 
1.67 % Xq introduced into the acceptance to measure- 
ments without converter. The converter increases the 
contribution from conversions and Dalitz decays by a 
fixed factor, which was determined precisely via GEANT 
simulations. Thus, the electron spectra from photonic 
and non-photonic sources can be deduced (Fig. EJ • The 
drawbacks of the converter method are the limitation in 
statistics of the converter run period and the fact that 
the photonic contribution is small at high pt- 

The e7 coincidence technique evaluates the correlation 
of electrons and photons via their invariant mass. Elec- 
trons from 7T° Dalitz decays or the conversion of one of 
the photons from 7r° — > 77 decays are correlated with a 
photon, in contrast to electrons from semileptonic heavy 
flavor decays. Comparing the measured e7 coincidence 
rate with the simulated rate for single 7r° events, allows 
to deduce Rnp as shown in Fig. once corrections for 
contributions from other photonic sources are applied. 

After subtracting the background cocktail from the in- 
clusive electron spectrum the invariant differential cross 
section of electrons from heavy flavor decays is shown 
in Fig. |3 compared with two theoretical predictions. 
A leading order (LO) PYTHIA calculation, tuned to 
existing charm and bottom hadroproduction measure- 
ments 23], is in reasonable agreement with the data for 
Pt < 1-5 GeV/c, but underestimates the cross section at 
higher pt- It is important to note that this calculation in- 
cludes a scale factor K = 3.5 to accomodate for neglected 
NLO contributions. A Fixed- Order plus Next-to- Leading- 
Log (FONLL) pQCD calculation [25j still leaves room 
for further contributions beyond the included NLO pro- 
cesses. The predicted contribution from bottom decays is 
irrelevant for the electron cross section at px < 3 GeV/c 
and becomes significant only for px > 4 GeV/c. 

The charm production cross section was derived 
from the integrated electron cross section for px > 
p Tlow = 0.6(0.8) GeV/c {da P e T - l ° w /dy = 4.78(2.15) ± 
0.78(0.46) (stat.) ± 1.74(0.68) (sys.) xl0~ 3 mb). Since in 
the low pt region, which dominates the total cross sec- 
tion, PYTHIA describes the measured spectrum reason- 
ably well, the total charm cross section was determined 
by extrapolating the properly scaled PYTHIA spectrum 
to px — GeV/c. First the PYTHIA spectra for elec- 
trons from charm and bottom decays were fit to the 
data for pt > 0.6 GeV/c, with only the normalizations 
as free parameters. The resulting midrapidity charm 
production cross section was determined as da c5 /dy = 
0.20±0.03(stat.)±0.11(sys.) mb, where the systematic er- 
ror is dominated by the uncertainty of the electron spec- 
trum itself (w56%), evaluated by refitting PYTHIA to 
the data at the minimum and maximum of the la system- 
atic error band. Additional uncertainties from the rela- 
tive ratios of different charmed hadron species and their 




'o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P T [GeV/c] 



FIG. 3: Invariant differential cross section of electrons 
from heavy flavor decays compared with PYTHIA LO (with 
K = 3.5) and FONLL pQCD calculations (upper panel). Er- 
ror bars (brackets) show statistical (systematic) uncertain- 
ties. For the FONLL calculation contributions from charm 
and bottom decays are shown separately. Ratio of data and 
FONLL calculation (lower panel) with experimental statis- 
tical (error bars) and systematic (brackets) uncertainties as 
well as the theoretical uncertainty (grey band). The solid line 
corresponds to the ratio of PYTHIA and FONLL. 



branching ratios into electrons («9%) and the variation 
of the PYTHIA spectral shape («11%) were added in 
quadrature. The rapidity integrated cross section was de- 
termined as a cS = 0.92±0.15(stat.)±0.54(sys.) mb, where 
various parton distribution functions (GRV98LO and 
MRST(c-g) [H in addition to the default CTEQ5L |3) 
were used for the extrapolation, with an associated extra 
systematic error of «6% [jj added in quadrature. 

Within errors the integrated charm cross section is 
compatible with data from Au + Au collisions 7\ (mini- 
mum bias value: 0.622 ± 0.057 ± 0.160 mb per NN colli- 
sion) and from d + Au collisions Q (1.3±0.2±0.4 mb) at 
the same y / sjv r Jv = 200 GeV. The FONLL cross section 
is smaller {a^ NLL = 0.256±°;^° mb) but it is still com- 
patible with the data. Our measurement does not allow 
to deduce a bottom cross section, which is predicted by 
FONLL as <J^ NLL = 1.87tg;$ fib. 

In conclusion, we have measured single electrons from 
heavy flavor decays in p + p collisions at = 200 GeV. 
These data provide a crucial benchmark for pQCD heavy 
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quark calculations. We observe that above pt ~ 2 GeV/c 
the electron spectrum is significantly harder than pre- 
dicted by a LO PYTHIA charm and bottom calculation. 
Contributions to the charm production cross section in 
excess of the considered FONLL calculation, e.g. from 
jet fragmentation, can not be excluded. The new data 
reported here provide an important baseline for the study 
of medium effects on heavy quark production at RHIC. 
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